Pseudogap and the Fermi surface in the t-J model 
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We calculate spectral functions within the t-J model as relevant to cuprates in the regime from 
low to optimum doping. On the basis of equations of motion for projected operators an effective 
spin-fermion coupling is derived. The self energy due to short-wavelength transverse spin fluctua- 
tions is shown to lead to a modified selfconsistent Born approximation, which can explain strong 
asymmetry between hole and electron quasiparticles. The coupling to long-wavelength longitudinal 
spin fluctuations governs the low-frequency behavior and results in a pseudogap behavior, which at 
low doping effectively truncates the Fermi surface. 
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Here we concentrate on some of the experimental facts 
revealing the nature of quasiparticles (QP) and pseu- 
dogap in undcrdoped cuprates. Several quantities, in 
particular the uniform susceptibility, the Hall constant, 
the specific heat, show the (large) pseudogap scale T* 
0, consistent with the angle resolved photoemission 
(ARPES) revealing a spectral function A(k, u>) with a 
hump at ~ 100 meV observed in Bi2Sr2CaCu202+,5 
(BSCCO) near the momentum k = (tt,0) §. QP dis- 
persing through the Fermi surface (FS) are resolved by 
ARPES in BSCCO only in parts of the large FS, in par- 
ticular along the nodal (0, 0)-(7T, tt) direction, indicating 
that the rest of the large FS are either fully or effectively 
gaped. Some aspects of the pseudogap have been found 
in the phenomenological spin-fermion models The 
renormalization group studies of the Hubbard model Q] 
also indicate on the breakdown of the standard Fermi liq- 
uid and on the truncation of the FS. That such features 
also emerge from prototype models of correlated elec- 
trons has been confirmed in numerical studies of spectral 
functions in the Hubbard || and in the t-J model HQ, 
which both show the appearance of the pseudogap at low 
doping. 

Our aim is to capture these features within an analyt- 
ical treatment of a single band model. The most difficult 
aspect in the latter is inherent strong coupling between 
mobile fermions and spin degrees, for which it is hard to 
find even a proper phenomenological model. In the fol- 
lowing we show that such an effective spin-fermion model 
can be derived via equations of motion (EQM) and di- 
viding the coupling into short and long-wavelength spin 
fluctuations. 

We study the planar t-J model 
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where we take into account possible longer range hop- 
ping, i.e., besides t%j = t for n.n. hopping also Uj = t' 
for n.n.n. neighbors on a square lattice. The latter ap- 
pears to be relevant in the study as shown later. We 
evaluate the single-particle propagator in this model ex- 
plicitly taking into account that fermionic operators are 
projected ones not allowing for double occupancy of sites, 
e.g., c| s = (1 — ni^ s )c\ s . We use EQM directly for pro- 
jected operators and represent them in variables rele- 
vant for a paramagnetic metallic state with (Sj) = and 
electron concentration (n.j) = c e = 1 — c^. EQM for £k s 
can be used to construct approximations for the electron 
propagator G(k,w) [fr|l0|], which can be represented as 



M~ Ck - S(k,w) : 
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where the renormalization a = (1 + c/J/2 is a conse- 
quence of a projected basis, and Ck = — 4r;it7k — 47^2^7^ 
is the 'free' propagation term emerging from the EQM, 
with 7k = (cos k x +cosfcj,)/2 and — cos k x cos k y . 
Here r/j = a + (SqSj) /a is determined by spin cor- 
relations. The central quantity for further considera- 
tion is the self energy E(k, uj) — ((Ck s ; C^ s )} l ^ r /a, where 
iCks = [ck s ,H] — CkCks, and only the 'irreducible' part 
of the correlation function should be taken into account 
in the evaluation of X. In finding an approximation for 
£ we assume that we are dealing with the paramagnet 
with pronounced AFM SRO with the dominant wave vec- 
tor Q = (tt, tt) and the AFM correlation length £ > 1 
with corresponding n = l/£. EQM naturally indicate 
on an effective coupling between fermions and spin de- 
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grees, however the role of short-range and longer-range 
spin fluctuations is quite different. 

In an undoped AFM system the spectral function of 
an added hole is quite well described within the self- 
consistent Born approximation (SCBA) pi, where the 
strong hole-magnon coupling induced by the hopping t- 
term leads to a broad background representing the inco- 
herent hopping and a quite narrow QP dispersion deter- 
mined by J. Our EQM formalism naturally reproduces 
SCBA in an undoped system and it is easy to generalize 
the equations for finite doping Ch > where we have also 
electron-like QP above the Fermi energy (ui > 0). In 2D 
the AFM long-range order is broken at T > and Ch > 0, 
still spin fluctuations are magnon-like with a dispersion 
w q for q > k and q > n where q = q — Q. The para- 
magnon contribution to the self energy can be written 
as 



E pm (k,a;) 
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[G~ (k - q, u + w q ) + G+(k + q, u - w q )] , 
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(1, -sign( 7q )) y/{2 J ± w q )/2w q and G ± 



refer to the Green's functions corresponding to electron 
(u > 0) and hole (w < 0) QP excitations, respectively. 

We are dealing with a paramagnet without the AFM 
long-range order, therefore it is essential to consider also 
the coupling to longitudinal spin fluctuations which for 
q < k appear to be quite uncoupled, therefore we express 
the longitudinal contribution as in Refs. §-0, 

^(k-q.u^x"^) 

UJ — 0J\ — U>2 

where x(Qj w ) is the dynamical spin susceptibility, 
A°(k,u) = -{ol/tz)Ivcx{u + fi - Ck - S pm ) _1 , 5(^1,^2) = 
[th(w 1 /2T) + cth(w 2 /2r)]/2and m kq = 2J 7q + ±(e°_ q + 
e°) with e° = -4t 7k - 4^0. 

In Sif only the part corresponding to irreducible dia- 
grams should enter, so there are restrictions on proper 
decoupling. We use in Eq. (^) the most appropriate 
and simplest approximation to insert the unrenormal- 
ized A°(k, w), i.e., the spectral function without a self- 
consistent consideration of Ey but with E pm fully taken 
into account. Such an approximation has been intro- 
duced in the theory of a pseudogap in CDW systems jl2) , 
used also in related works analyzing the role spin fluctu- 
ations ||l3|| , fy|] , and recently more extensively examined 
in Ref. IJ. 

So far we do not have a corresponding theory for the 
spin response at Ch > and T > 0, so x(k, oj) is assumed 
as a phenomenological input, bound by the sum rule 

rcth(^)x"(q, W )d«=|(l-c h ). 
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The response should qualitatively correspond to a para- 
magnet close to the AFM instability, so it is assumed of 
the form 



X"(q,w) oc 
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where 4>(uj, T) oc u> would be appropriate for a nearly 
AFM Fermi liquid [g,[l4) or an undoped AFM in 2D. 
The selfconsistent set of equations for G is closed with 
E = E pm + Eif. For given chemical potential \x the 
FS emerges as a solution determined by the relation 
Ck F + E'(ki?,0) = /1. We should note that at given /z, 
electron concentration c e as calculated from the density 
of states Af(w) = (1/N) £) k A(k, lu) integrated for uj < 0, 
does not in general coincide with the one evaluated from 
the FS volume, c e = Vfs/Vq. Nevertheless, apart from 
the fact that within the t-J model validity of the Lut- 
tinger theorem is anyhow under question |lq| , in the 
regimes of large FS both quantities appear to be quite 
close. The position of the FS is mainly determined by 
and E pm , while in this respect Ejf is less crucial. We 
choose further on parameters J = 0.3t, t' — — Q.2t and 
k = Jch while 771 and 772 are determined as a function 
of Ch from model calculations 0|. We use N — 32 x 32 
points in the Brillouin zone and broadening e/t = 0.02. 

In Fig. 1 we present hole concentration Ch vs the chem- 
ical potential fi/t as obtained from Eq. 0. We solve 
the equation by iteration. With labels (a), (b) and (c) 
are indicated special cases presented in next figures. At 
Ch ~ 0.12 we observe in the equations an instability sig- 
naled by oscillatory behavior instead of convergence and 
a unique solution can not be obtained in the region in- 
dicated by the dashed line. However, at lower doping, 
Ch < 0.05, the solution converges again and a typical 
result is indicated with (d). The region of instability co- 
incides with the transition from the large to a small FS, 
as presented below. 

In Fig. 2 are presented spectral functions A(k, uj) along 
the principal directions in the Brillouin zone for cases 
(a), (c) and (d) from Fig. 1. It is evident that E pm 
leads to strong damping of hole QP and quite incoherent 
momentum-independent spectrum A(k, ui) for lu <C —J 
which qualitatively reproduces ARPES and numerical re- 
sults . Electron QP (at u> > 0) are in general very dif- 
ferent, i.e., with much weaker damping arising only from 
E pm . At low doping Ch < 0.05 we find the regime of small 
(pocket-like) FS. 

The shape of the FS is most clearly presented with con- 
tour plot of the electron momentum distribution function 
n(k) = a -1 f_ A(k,Lu)duj. Results for a characteristic 
development of the FS with are shown in Fig. 3. At 
Ch < c*h ~ 0-05 solutions are consistent with a small 
pocket-like FS, (d), whereby this behavior is enhanced 
by t' < as realized in other model studies p^j. On 
increasing doping the FS rather abruptly changes from 
a small into a large one as suggested from the results of 
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SCBA jL!|. The smallness of c* h has the origin in quite 
weak dispersion dominated by J and t' at Ch — > which 
is overshadowed by much larger £k at moderate doping, 
where the FS is large and its shape is controlled by t'/t. 
Figures (a), (b) and (c) correspond to higher doping with 
common large FS topology. However, in the intermediate 
doping regime, (b) and (c), the pseudogap is pronounced 
at momenta around (tt, 0) point. The gap is more pro- 
nounced in (c) because of longer AFM correlation length 
£ (smaller k). In Fig. 4 the development of the PG is 
presented for Ch = 0.21 [case (b)]. Here k = 0.4 and the 
PG is not fully developed yet. With increasing £ - for 
this doping somewhat unrealistic values - the gap opens. 
At lower doping, case (c) in Fig. 3 the gap is opened what 
also shows up the form of the FS which tends to avoid 
points (ir, 0), (0, 7r). 

In conclusion, we have presented a theory for the spec- 
tral functions within the t-J model where the double- 
occupancy constraint is taken explicitly into account and 
used to derive an effective spin-fermion coupling. The 
coupling to transverse AFM paramagnons is strong, nev- 
ertheless it can be well treated within a generalized 
SCBA. On the other hand, the coupling to longitudi- 
nal AFM fluctuations, liikq, is moderate near FS for low 
doping and leads to a pseudogap. The latter is however 
not in contradiction with the existence of a large FS, and 
should show up in integrated photoemission and ARPES 
results as well as in the uniform susceptibility and in the 
specific heat. 



44, 317 (1991); A. Ramsak and P. Prelovsek, Phys. Rev. 

B 42, 10415 (1990). 
[12] PA. Lee, T.M. Rice, and P.W. Anderson, Phys. Rev. 

Lett. 31, 462 (1973). 
[13] A. P. Kampf and JR. Schrieffer, Phys. Rev. B 41, 6399 

(1990). 

[14] A.V. Chubukov and D.K. Morr, Phys. Rep. 288, 355 
(1997). 

[15] A.J. Millis and H. Monien, Phys. Rev. B 61, 12496 
(2000). 

[16] W.O. Putikka, M.U. Luchini, and R.R.P. Singh, Phys. 

Rev. Lett. 81, 2966 (1998). 
[17] J. Bonca, P. Prelovsek, and I. Sega, Europhys. Lett. 10, 

87 (1989). 

[18] T. Tohyama and S. Maekawa, J. Phys. Soc. Jpn. 59, 1760 
(1990). 

[19] A. Ramsak, I. Sega, and P. Prelovsek, Phys. Rev. B 61, 
4389 (2000); A. Ramsak, P. Prelovsek and I. Sega, in 
Proceedings of the NATO ARW on Open Problems in 
Strongly Correlated Electron Systems, edited by J. Bonca 
et al. (Kluwer, Dordrecht, 2001). 



FIGURES 



[1] for a review see e.g. M. Imada, A. Fujimori, and Y. 

Tokura, Rev. Mod. Phys. 70, 1039 (1998). 
[2] D.S. Marshall et al., Phys. Rev. Lett. 76, 4841 (1996); 

H. Ding et at, Nature 382, 51 (1996). 
[3] J. Schmalian, D. Pines, and B. Stojkovic, Phys. Rev. 

Lett. 80, 3839 (1998). 
[4] D. Zanchi and H.J. Schulz, Europhys. Lett. 44, 235 

(1997); N. Furukawa, T.M. Rice, and M. Salmhofer, 

Phys. Rev. Lett. 81, 3195 (1998) 
[5] R. Preuss., W. Hanke, C. Grober, and H.G. Evertz, Phys. 

Rev. Lett. 79, 1122 (1997). 
[6] J. Jaklic and P. Prelovsek, Phys. Rev. B 55, R7307 

(1997); Prelovsek, J. Jaklic, and K. Bedell, Phys. Rev. 

B 60, 40 (1999). 
[7] for a review see J. Jaklic and P. Prelovsek, Adv. Phys. 

49, 1 (2000). 
[8] P. Prelovsek, Z. Phys. B 103, 363 (1997). 
[9] N.M. Plakida and V.S. Oudovenko, Phys. Rev. B 59, 

11949 (1999). 

[10] P. Prelovsek and A. Ramsak, Phys. Rev. B, 180506(R) 
(2001). 

[11] C. L. Kane, P. A. Lee, and N. Read, Phys. Rev. B 39, 
6880 (1989); G. Martinez and P. Horsch, Phys. Rev. B 



0.3 



0.2 



0.1 



0.0 



— i — r 


/ 


—i — r 


— i 1 — | 1 — i - 




I i i i i I i i i i 












J=03t 




(a) 






f=-o.2f ; 








(b) 
























(c) 


^\ r | - 



0.0 0.5 1.0 1.5 2.0 2.5 



\llt 

FIG. 1. Hole concentration Ch as a function of the chemical 
potential fi/t for J/t = 0.3, t'/t = -0.2 and N = 32 x 32 sites 
in the Brillouin zone. Note selected cases presented further 
on: (a) c h = 0.28, (b) c h = 0.21, (c) c h = 0.16 and (d) 
c h = 0.04. 
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FIG. 2. (color) Spectral functions A(k.,ui) along principal lines in the Brillouin zone. For convenience of presentation the 
functions are clipped at A(\c,u) = 1.0. Labels (a) - (d) correspond to selected cases in Fig. 1. 
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FIG. 3. (a,b,c,d) (color) Electron momentum distribution function n(k) for various as selected in Fig. 1. 
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